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Abstract 

We establish a variant of the symmetric quasi linear first order 
system given by H. Priedrich for the evolution equations of gravi- 
tating fluid bodies in General Relativity which can be important to 
solve realistic problems. Our version has the advantage of introducing 
only physical characteristics. We state explicitly the conditions un- 
der which the system is hyperbolic and admits a well posed Cauchy 
problem. 

Dedicated to Andrzej Granas. 



1 Introduction. 

Cattaneo [C] and Ferrarese [Fl], have stressed long ago that the fundamen- 
tal physical data are those denned on time lines instead of those defined on 
spacelike hypersurfaces. In the case of a fluid, timelike world lines have a 
material reality given by the flow, the unknowns associated to them corre- 
spond to a lagrangian picture, in contrast with the usual eulerian one. In 
many situations, particularly in cosmology, the lagrangian description seems 
more fundamental. However, one cannot be too confident because General 
Relativity holds many surprises. For example, consider a fluid in internal 
equilibrium and in equilibrium with a stationary black hole. Under some 
conditions (cf. [BY]) the thermodynamic temperature is found by certain 
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eulerian observers whose world lines have vanishing vorticity. This corre- 
sponds to our "well posedness" result given in section 11.4. We believe that 
eulerian and lagrangian descriptions are both usefull. 

Ferrarese [F2], [F3] has given a complete 3+1 decomposition of the con- 
nection and Riemann tensor in a frame whose time axis is tangent to the time 
lines and the space axes orthogonal to it. Using this decomposition he has 
written the Einstein- dust equations as a (infinite dimensional) dynamical 
system. Recently H. Friedrich [HF] has used a lagangian description of the 
flow together with a tetrad formalism to write the Einstein equations coupled 
with a perfect fluid as a symmetric first order system, formally hyperbolic. 
His treatment includes dust as a particular case. 

The Einstein- dust system, as well as the Einstein- Euler system for a 
perfect fluid, have been proved long ago [C.B1] to form a well posed hyper- 
bolic system in the sense of Leray. The relativistic Euler system has been 
put in first order symmetric hyperbolic form (FOSH) first by K.O. Friedrichs 
himself [KOF] , then by general methods by Anile, Boillat and Ruggeri (see 
references in [A]). Rendall [R] has obtained, for well chosen equations of 
state, FOSH systems for the Euler equations which extend to vacuum, but 
do not yet permit, as he points out himself, the solution of realistic problems 
of motion of compact fluid bodies in otherwise empty space. The system 
obtained by H. Friedrich can be important to solve such realistic problems. 

In this article we establish a variant of Friedrich's system for a lagrangian 
and tetrad description of the flow and the spacetime geometry. Contrary to 
Friedrich we will not work with the Weyl tensor but with the Riemann tensor, 
as in [AC. BY]. It avoids the introduction of unphysical characteristics and 
seems to us more natural in this non vacuum case as well as somewhat sim- 
pler. The systems we obtain for vacuum, dust or perfect fluid are symmetric 
and hyperbolic in the same sense as Friedrich's system. They are hyperbolic 
in K.O. Friedrichs's sense, which leads to existence theorems for solutions 
in Sobolev spaces, only if the time lines admit global spacelike sections, and 
the initial data are given on such a section. These existence theorems are 
needed even if one wants to consider the evolution equations as a well posed 
dynamical system. 
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2 Metric and coframe. 

We write the metric in an orthonormal frame, i.e., 

3 

9 = -(0°) 2 + J2^ 2 - (2- 1 ) 

i=i 

We choose the time axis to be tangent to the time lines, i.e. the cobasis 
9 is such that Q % does not contain dx°. We set 

6 l = a)dx j , 6° = Udx° + bidx\ (2.2) 

We will call such a frame a C.F. (Cattaneo- Ferrarese) frame. We denote by 
A\ the matrix inverse of a\. It holds that 

dx i = a^, dx° = u-\e° - A)b t e r ) . 

The Pfaff derivatives d a in the C.F. frame are linked to the partial deriva- 
tives d/dx a by the relations 

d d 



dx° ' dx 

3 Bianchi equations 

Whatever the coframe, the components of the Riemann tensor satisfy the 
identities 

^ a R/3y,\fi + VpRja^n + V 7 -R a/ 3 j A A t = 0, (3.1) 

hence, if the Ricci tensor R a/ 3 satisfies the Einstein equations 

R a(B = p a p, (3.2) 

it holds that 
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Equations 3.1 and 3.3 imply the following ones (cf. an analogous system in 
[AC. BY]) 

Vo-R/ifc,oj + Vfc-Ro/i.oj ~~ Vft-Rofc,oj = 0, (3.4) 
and, in using a symmetry of the Riemann tensor, 

+ = Ji,oj = V oPji - V jPoi . (3.5) 

The equations 3.4 and 3.5 are for each given pair (Oj) a first order system 
for the components Rhk,oj and Roh,oj- The principal operator is a symmetric 
6 by 6 matrix: 





do 








d 2 







\ 







<% 








9 3 


-<9 2 












do 


-9 3 





& 






d 2 





-d 3 


9 












-di 










<9 







V 





-9 2 


9! 








do 


) 



Analogous equations and result hold with the pair (Oj) replaced by (Im), 
I < m. The system finally obtained has a principal matrix consisting of 6 
identical 6 by 6 blocks around the diagonal, it is symmetric. 

However, the system obtained for the components of the Riemann tensor 
in the coframe with time axis 9° and space axis 6\ with principal operator 
the matrix M a d a , cannot be said to be hyperbolic in the usual sense: the 
principal matrix M° is the unit matrix hence positive definite, but the oper- 
ators d a are not partial derivatives. We say that the system is a quasi FOSH 
(First Order Symmetric Hyperbolic) system (Cf. section 11). 

Remark 1 One can associate with the Riemann tensor, as in the usual case 
of the 3+1 splitting, two pairs of "electric" and "magnetic" 2- tensors whose 
non zero components in the frame are: 

rp — p 7~) = }_c c r>hk,lm 

Hilj — ±\,QiQj, Uij — ^ihk^jlm 1 ^ 

TT — }_ c r>hk r> — }_ c r>:::::hk 

12 ij — 2 ihk Oj' ji — ^ihk^Oj ' 

where is the totally antisymmetric Kronecker tensor. These definitions 
imply here 

Rhkfij = e \;hkHiji Rhk,lm = e \:hk E \:l-mP ij ■, Roj,hk = £ l.hk^ji- 

The previous system is a quasi FOSH system for these pairs of tensors. 
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Remark 2 The system 3.1, 3.2 contains also the equations, which we do not 
use for evolution, 

^hRij,\n + VjRhi,\[i + ViRjh,\n = (3.6) 

and 

V h R h 0M = V aPm) - V mPao . (3.7) 

These equations are not usual constraints on a submanifold t = constant 
because di contains the transversal derivative d/dt. We call them Bianchi 
quasi- constraints . 

4 Coframe structure coefficients. 

The structure coefficients c of the coframe are defined by 

dd a = -\c a Pl e p A e\ 

We have 

de° = diUe 1 a dx° + d a be a a dx l 

that is, 

d9° = U- l d t U9' t A (9° - A h jb h 9 3 ) + A h j d a b h 9 a A 9 j , 

therefore 

cl = U-%U - Ajd b, = -c% 
and, with f {ij] = fa - f jh 

4 = U- l b h A^U - A h b d i} b h 

that is 

<*=UAfa< l (U- 1 b h ) 

also 

d9 l = A{(d h a)9 h A9 k + d o a)9 A 9 k ) 



5 



hence 



4o = Aidoa) 



4k = A l h dk]a) 



Remark 3 If the time lines are not hyper surface orthogonal ( i. e. if hi 7^ 0) 



the coefficients cf- are different from the structure coefficients of the space 



frame 9 l . 

5 Connection. 

The connection is such that 



LUfi,^, — U^o — Cc 



and 

W 7/3,A + w 7 A,/3 = 

with (77 is the Minkowski metric) 
We set 

ol 

C 7 /3,A = VaX C yPi 

and we have (Cf. C.B-DeWitt I p.308) 

We find that, as foreseen from antisymmetry 
We set 

Yi = o;oo,i, 
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and we have 

Yi = u 00 ,i = uJq { = u 00 = -coi,o = c oi = U^diU - A\d$j. 
We know that uj ] Qi = u 0i j is antisymmetric in i and j. We set 

u 0hJ = hi = \{A h 3 d^ h - A%al + 4^M- 

Remark 4 Let e( a ) = <9 a fre the frame dual to 9 a , i.e. with components S^ a y 
Then 

in particular uoij is the projection on of the derivative of e^ in the 
direction of e( ). We have fij = if the frame is Fermi transported along the 
time line. 

The connection coefficient ujio t j is the sum of a term symmetric in i and 
j and an antisymmetric one, we have 

X l3 = uj i0 , 3 = col, = ^{Ap ai + A%a\ + A^b h }. 

The antisymmetric term vanishes if the time lines are hypersurface orthogo- 
nal (h = 0). 

The coefficients lu^ are also linear expressions in terms of the first deriva- 
tives of the frame coefficients, they are identical to the connection constructed 
with the a\ at fixed x° if b, = 0. 

6 Frame evolution. 

The relations between the connection coefficients and the frame coeffi- 
cients a\ give equations for the dragging of these coefficients along the time 
lines when the connection is known. Indeed we have found 

d a) = aJ4 = 2a}(a4 - ^ ) = 2a*tfj - X*) (6.1) 
dob, = a£(-Y h + U-'dnU (6.2) 
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No equation gives the evolution of U. It can be considered as a gauge 
variable fixing the time parameter. H. Friedrich chooses the coordinate x° to 
be the proper time along the given time lines, i.e. U = 1 

The quantity /j is also a gauge variable fixing the evolution of the space 
frame. If we choose it to be Fermi transported like Friedrich, then /j = 0. 

7 Curvature. 

7.1 Evolution of connection. 

Using the general formula (cf. C.B- DeWitt I p. 306, with opposite sign 
convention) 

Rx^p = dxufo - d^y - - u^) + - u«y xp (7.1) 

R'uh'.j = d^\ 3 - V h f] -Y h f]-Uh- XK k W kj + Y l X hj - YjX^ (7.2) 

where V is the pseudo covariant derivative constructed with <9j and 
(Cattaneo- Ferrarese transversal derivative). In the Fermi gauge Rq^j re- 
duces to 

Roh^ = douij + XMj + W ~ Y,X h \ (7.3) 
The principal operator in these equations is the dragging along the time lines 

of 4. 

We have: 

Rmio = V;,V, - OoXjf + ) ) ), - XjXf + tfX? - tfxj. (7.4) 

We deduce from this identity, in the Fermi gauge, 

R 00 = W l Y i - d Xi + Y% - X- h j X- h . (7.5) 

In the next sections we will write other equations to achieve the determi- 
nation of both Y and X. 
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7.2 Quasi-constraints. 

The other components of the Riemann tensor are given in the C-F frame by: 

R hk!.j = R hk':..j - f}( x hk - x kh ) + x-j.x jh - x jk x l h , (7.6) 

where R'^.j denotes the expression formally constructed as a Riemann tensor 
with the coefficients a;?,-. 

R k T°, = d k u° hj - d h ul - u° pj (^ kh - u p hk ) + uly hj - uly kp (7.7) 

that is ,with previous notations, 

R'kh.'°.j = ^kXhj — VhX k j — Yj(X kh — X hk ). (7.8) 

These equations do not enter in the evolution system for the connection 
that we are considering, they do not contain the derivative <9o, we call them 
connection quasi- constraints. 

Remark 5 Modulo the expression of the connection in terms of frame coef- 
ficients one has the well known symmetry: 

Rkhflj = Roj,kh 

One deduces from the identity 7.8 

R h0 = VjX^ - V h X] - Yi(X jh - X hj ) 

8 Vacuum case. 

In the vacuum case there are no a-priori given time lines. We can give arbi- 
trarily on the spacetime the projection Yi of V e(0) e(o) on e^y The quantities 
U and f\ being also arbitrarily given the Bianchi equations together with 
6-1, 6-2; 7-3; 7-5 constitute a quasi FOSH system for the unknowns Eij, Dij, 
Hij, Bij, aj, bi, u^, Xij. Its solution determines the spacetime metric. 
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9 Perfect fluid. 

9.1 Fluid equations. 

The stress energy tensor of a perfect fluid is 

T a/3 = (/i + p)u a up + pg aP . 

Then 

One supposes that the matter energy density fi is a given function of the 
pressure p and entropy S; this last function is conserved along the flow lines 

u a d a S = 0. (9.1) 

The Euler equations of the fluid express the generalized conservation law 

W a T aP = 0, 

and are equivalent to the equations 

(H + p)u a V a uP + (u a v? + g a ^)d aP = 0, with u a u a = -1 

and 

(fx + p)V a u a + u a d a fi = 
In our coframe they reduce to: 

d S = (9.2) 

(ti+p)Y i + d i p = 0, 1^ = ^00 (9.3) 

d of i+(fi + p)Xl = 0. (9.4) 
Using the index F of the fluid defined by 

dp 



F(p, S) 



fi(p, S)+p' 
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the equation 9.3 reads 

Yi = -diF, (9.5) 

while 9.4 gives , modulo 9.2, 

fi' p d F + XI = 0. (9.6) 
The commutation relation between Pfaff derivatives 

(d d l - d t d )F = c^d a F 

implies therefore 

^[doYi + Yd F + (// - Xi)djF] - d^doF - d,X h h = 0. (9.7) 

The use of 9.5 and 9.6 replaces d a F by functions of Y,X and p,S. The 
derivatives dip! p are functions of Y,p, S and dkS, since 

difi'p = n" p 2dip + n" pS diS. 

We introduce Sk = d^S as new unknowns. They are dragged along the flow 
lines by the following equation deduced from 9.2: 

daS^^Sj^ifi-X})^. (9.8) 

Following H. Friedrich we replace in 9.7 diX% by its expression deduced 
from the equation 

R i0 = V h X-- h - ViXjJ - Y h (X hi - X th ) = 0, 

and we obtain, changing names of indices 

fi' p d Y h - V 3 X{ - Y\X hl - X lh ) + fi' p [Y h d F + 

(f h - X{)d 3 F] + d h »' p d F = 0. (9.9) 

In 7.4 we replace VhXi by V^Y^ + c^Yq, with 

Y = -d F = -(fjQ-'Xi c ° hi = X ih -X hi 

and we obtain 

doXjl - vV fc - c° hi Y - )),) ' + i.\ h l - ti)Xf + jyX h > = -R^ . (9.10) 

The principal operator on the unknowns Y and X in the equations 9.9, 
9.10 is diagonal by blocks. Each block, corresponding to a given index h, is 
symmetric, it reads: 
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If fi' v > the matrix M° is positive definite in the C.F. frame. The system 

j 
h- 



9.9, 9.10 is a quasi FOSH system for the pairs Yh,Xf 



Remark 6 The characteristic determinant associated with the system 9.9, 
9.10, obtained by replacing d a with a covariant vector C, is 

KM- E ^ 3 

i=l,2,3 

The roots of — Xli=i2 3^i = correspond to sound waves. Their speed 
is at most 1 (speed of light) if and only if fi' p > 1. 

The full system of fluid equations is quasi diagonal by blocks. The ma- 
trices M a d a corresponding to 9.2, 9.4, 9.8 reduce to <9 . 

9.2 Sources of the Bianchi equations. 

In our frame the source tensor p a/3 reduces to 

Poo = \{» + 3 P)> Poi = 0, Pa = \^- P) 6 v 

We have seen that d a p and d a /i are smooth functions of p, S,Si,Y and X. 
The same property holds for J it0 j and J^hk- 



9.3 Conclusion. 

Assembling the results of the previous sections we find the following theorem. 

Theorem 7 The EEF ( Einstein- Euler-Friedrich) system for a gravitating 
perfect fluid, with flow lines taken as time lines, U and /j = uj Q3 - given arbi- 
trarily, is a quasi FOSH system for the Riemann curvature tensor, the frame 
and connection coefficients, the density of matter, the entropy and its space 
derivatives. 

As remarked in section 3 the system is not a usual FOSH system: the 
time lines are by choice timelike but the hypersurfaces x° = constant are 
not necessarily spacelike for the characteristic cone. 
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10 Case of dust. 



10.1 Matter equations. 

The stress energy tensor of a dust source is 

The flow lines, tangent to the unit vector u a , are then geodesies, it holds that 

u a V a vP = 0, (10.1) 
while the conservation of matter reads: 

u a d a [i + nV a u a = 0. (10.2) 
We take the flow lines as time lines, i.e. u a = 5%. Then 10.2 reads 

d /i + nX\ = 0, (10.3) 

while 10-1 gives 

c4o =Yi = 0. (10.4) 

Using this equation we see that, given arbitrarily on the spacetime the gauge 
variables a[j the equations (6-1), (6-2) reduce to a quasi- diagonal system with 
principal operator d for the frame coefficients when the connection (which 
appears undifferentiated) is known . 

10.2 Bianchi equations. 

For a dust stress energy tensor it holds that 

Pa/3 = Ku a Up + -g a p), 

hence in the chosen frame 

POO = 77P; Poi = °> Pij = 
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A straightforward computation gives then, choosing the Fermi transported 
frame for simplicity 

Ji,oj = VoPjj - Vjp 0i = fj,(--5ijX% + Xji) 

Ji,hj = ^hPji - VjPhi = ^{dhpSji - djfi5 hi + fi(u l hj - u) h )} 

There appears a difficulty which is the presence of the space derivatives 
dh\x in Ji t hj- We get rid of this problem, inspired (though with a different 
choice) by H. Friedrich's treatment, by replacing the unknown R'^.j by an- 
other unknown with the same symmetries, namely 

Fij t hk. = Rij,hk + ^ikRjh — $ihRjk — SjkRih + $ihRik + Dij,hkR 

with 

Dij,hk = -^{o~ik$jh — &jk&ih) 

Using this unknown, the contracted Bianchi identities become 

V R° ::i , hj + Vk !■%.,„ = 0, (10.5) 
while the original ones become 

Vo-Fjj./ife + VjRoi,hk — Vj-Roj>fc = 

Vo{5ikRjh — SihRjk — SjkRih + SihRik + Xij hk R} (10.6) 

One obtains a quasi FOSH system for Roi,hk and F^hk by replacing on the 
right hand side of this equation R a p by p aj3 : now there appears only the 
time derivative <9 /i which, satisfying 10-3, can be eliminated in favor of 
undifferentiated terms. 

Note that Rij^k does not appear in other equations than the Bianchi 
equations. 

10.3 Conclusion. 

We have proved the following theorem 

Theorem 8 The EDFF (Einstein-Dust-Ferrarese-Friedrich) system, , with 
flow lines taken as time lines, U and /j = uo l - given arbitrarily, is a quasi 
FOSH system for the Riemann curvature tensor, the frame and connection 
coefficients, the density of matter and its space derivatives. 



14 



11 Cauchy problem. 



11.1 Hyperbolicity. 

The fact that for the systems that we have obtained the hypersurfaces x° = 
constant are not necessarily spacelike for the characteristic cones poses a 
difficulty for the Cauchy problem, since energy estimates used in proving 
existence theorems rely on this property. 
We have, using standard definitions 

Proposition 9 A quasi FOSH system, with principal operator the matrix 
M a d a , is a FOSH system with x° as a time variable if the matrix of coeffi- 
cients of d/dx° is positive definite, namely if 

M° = M a A° a = U'\M° - A)biM j ) 

is positive definite. 
11.1.1 Case of dust 

Lemma 10 The EDFF (Einstein-Dust-Ferrarese-Friedrich) system is a FOSH 
system relative to x° = cons tan t slices as long as the quadratic form 

9 jh = a )<-bjbh (H.l) 

i=l,2,3 

is positive definite and U > 0. 
Proof. 

The matrix M^ bian ^ of coefficients of d/dx° corresponding to the Bianchi 
equations is after multiplication by U, setting 









B t = 
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Its eigenvalues are, each with multiplicity 2: 

1, 1 + ^{Bl + Bl + B\\ 1 - yJ{Bl + Bi + B\) 

They are positive, because the given condition implies B\ + B\ + B\ < 1. 
The other matrices M?, and M? tt , are unit matrices. 

11.1.2 Fluid case. 

Lemma 11 The EEF ( Einstein- Euler-Friedrich) system is a FOSH system 
relatively to x° = cons 1 slices as long as the quadratic form gjh given in 11.1 
is positive definite, U > and fj/ > 1. 

Proof. The principal operator is symmetric and composed of blocks 
around the diagonal. 

The matrix M? U s is the same as in the case of dust. 

The matrix M° on corresponding to the connection evolution is, after mul- 
tiplication by U 



( 


/4> 


-B 1 


-B 2 


-B 3 


\ 




-B l 


1 












-B 2 





1 
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-B 3 
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/ 



Its eigenvalues are 1, with multiplicity two, and: 



1 

2 



;{l + /4±2^(^-l) 2 + ^ 2 )} (11-2) 
These eigenvalues are positive under the hypothesis 11-1 if ji' p > 1. 



Remark 12 The condition ji' p > 1 expresses that the sound speed is at most 
equal to the speed of light. 

11.2 Cauchy data. 

The Cauchy problem is less natural in the C.F. formulation than in the usual 
3+1 decomposition. 
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11.2.1 Initial data. 



An initial data set on a coordinate patch of a 3 dimensional submanifold is 
composed of the following elements: 

1. A field of coframes a*- 7 and a field of covariant vectors aS 
The quadratic form on M with coefficients 

9jh= ^2 affi-bjbh (11.3) 

i=l,2,3 

is supposed to be positive definite. The submanifold M = M x {0} of M x R 
is then spacelike for any lorentzian metric on M x R which reduces on M to 

-(Udx° + bidx*) 2 + (^) 2 . U > 

i=l,2,3 

2. Quantities u)^-, A/,Fj. 

3. Pairs of 2-tensors E % \ H^, D^, By. 

4. a (case of dust) A scalar function p, > 0. 

4.b (perfect fluid case). Two scalar functions p > and 5 > 0, and 
equation of state /i(p, S) such that /i' p > 1 for all p in a neighbourhood of 

(p,S). 

11.2.2 Constraints. 

The initial data on M do not determine on M the derivatives <9j = di—Albjdo. 
To satisfy the quasi-constraints on M we must first deduce from the given 
data and the considered evolution system the values on M of the derivatives 
with respect to x° of the relevant unknowns. The quasiconstraints give then 
relations on the initial data which we call constraints and suppose to be 
satisfied on M. 

11.3 Existence theorem for the reduced system. 

Local existence theorems for a solution of the Cauchy problem are a direct 
consequence of known results for quasilinear FOSH systems. We enunci- 
ate the result in the perfect fluid case. The spaces Hg doc (M) are the usual 
uniformly local Sobolev spaces on M. 
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Theorem 13 Let there be given an initial data set in H^" loc (M), s > § + 1. 
Then for any choice on M x I of the gauge variables in C°(I, Hg doc (M)),with 
I an interval of R containing 0, there exists a neighbourhood I' ofO in R such 
that the EEF equations admit a solution in C°(I', iJ s M - /oc (M))nC 1 ( // > H^° C (M)) 

To prove that the solution satisfies, modulo the initial constraints, the 
original Einstein- Euler equations would be a messy task, which we will not 
endeavour. The uniqueness of the obtained solution gives an indirect proof, 
modulo the existence theorem already known [C.B1], [C.B2]. 

11.4 Irrotational flows. 

We recall the equation 6-2: 

d b t = a1(-Y h + d h logU), 
with, in the case of a perfect fluid, 

Y h = -d h F. 

We choose then (gauge choice, the geometrical result is independent of it) 

U = exp(-F). 

The equation reduces then to 

d bi = 0. (11.5) 

Therefore it holds that b{ = througout the motion if it is so initially. It is 
the well known property of conservation of zero vorticity (irrotationality) of 
a perfect fluid flow. 

For an irrotational flow it holds that is a symmetric 2- tensor, equal 
up to the sign to the extrinsic curvature of the space slices x° = constant 
(components in the frame 9 l ), and = U~ 1 diU, while = lv^ are the 
connection coefficients of the space metric g = X]«=i2 3(^) 2 ' ® % = a h ( ^ xh - 
The pseudo covariant derivatives Vj are covariant derivatives V« in the space 
metric. 

For an irrotational flow the EFF system is a FOSH system. 
The same result holds in the case of dust the appropriate gauge choice is 
then U = constant. 



(11.4) 
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Remark 14 The relations between frame and connection, as well as connec- 
tion and curvature, for hyper surf ace orthogonal time lines are a special case 
of those found in the 3+1 decomposition [C.B3J. In particular, in the present 
notations, the Einstein constraints take the familiar form: 

R l0 = VjX{ - diX] = (11.6) 

S 00 = i(i2oo + PR*) = \{R ~ Xjxf + (XI) 2 ) = Poo = + 3p). (11.7) 
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